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Spatio-temporal dynamics of a deterministic three-level cellular automaton (TLCA) of Zykov- 
Mikhailov type (Sov. Phys. - Doklady, 1986, Vol.31, No.l, P.51) is studied numerically. Evolution 
of spatial structures is investigated both for the original Zykov-Mikhailov model (which is applica- 
ble to e. g. Belousov-Zhabotinskii chemical reactions) and for proposed by us TLCA, which is a 
generalization of Zykov-Mikhailov model for the case of two-channel diffusion. Such the TLCA is a 
minimal model for an excitable medium of microwave phonon laser, called phaser (D. N. Makovetskii, 
Tech. Phys., 2004, Vol.49, No. 2, P.224; jcond-ma t/0402640 ). The most interesting observed mani- 
festations of TLCA dynamics are as follows: (a) spatio-temporal transient chaos in form of highly 
bottlenecked collective evolution of excitations by rotating spiral waves (RS Ws) with variable topo- 
logical charges; (b) competition of left-handed and right-handed RSWs with unexpected features, 
including self-induced alteration of integral effective topological charge; (c) transient chimera states, 
i. e. coexistence of regular and chaotic domains in TLCA patterns; (d) branching of TLCA states 
with different symmetry which may lead to full restoring of symmetry of imperfect starting pattern. 
Phenomena (a) and (c) are directly related to phaser dynamics features observed earlier in real 
experiments at liquid helium temperatures on corundum crystals doped by iron-group ions. 
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I. INTRODUCTION 

During last years vortices and, in particular, rotating 
spiral waves (RSWs) are widely studied in excitable me- 
dia [S HI 12, lasers IH El Hi IH 113 , Bose-Einstein 
condensates |2a|, saturable nonlinear media (£3, ferro- 
magnetics |3ll l59j|. photorefractive optical lattices |96| . 
open microwave billiards |3 5 j. th ree- level spatial Lotka- 
Volterra's type models 

[lOflUoiJ etc. 

Evolution of such complex spatio-temporal structures, 
as patterns of RSWs, is now a subject of extensive com- 
puter modeling, especially for intrinsically unstable sys- 
tems possessing hypersensitivity to initial and/or border 
conditions, to imperfections of active medium etc. 

Choice of appropriate model is the crucial point for 
computer studying of these complex systems. Even fully 
deterministic model of continuous unstable system may 
demonstrate unpredictable behaviour due to exponential 
increasing of small deviation of initial conditions from 
their precise values |4lj . This unpredictability cannot 
be eliminated because of finite precision of initial data, 
available for a digital computer, even if no round-off er- 
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rors exist in certain "perfect" digital machine. Generally, 
instead of computing of the evolution of a continuous un- 
stable system (in the framework of given model), a digital 
machine may compute an evolution of some composite 
system, which includes the computer itself. 

Another serious trouble is the large dimension prob- 
lem, which is typical for modeling of many-particle sys- 
tems. Using of appropriate averaging is a good solution 
for systems near equilibrium, but for strongly nonequlib- 
rium systems it often is not so. 

A possible way to bypass these troubles consists in us- 
ing of part ially or fully discrete analogs of continuous 
models [jlGlMlIlllllIllIlS Con- 
struction of such models is very nontrivial task, but ade- 
quate discrete model is much more suitable for numerical 
investigation of unstable and/or multiparticle systems, 
than, say, a model based on partial differential equations. 

In this work, we study evolution of RSWs and some 
other spatio-temporal structures in discrete active (ex- 
citable) system consisting of locally interacting three- 
level particles (units) with finite times of relaxation. This 
active system is described by Zykov-Mikhailov (ZM) mo- 
del, which was primarily introduced in 1986 for chemical 
excitable systems [9^ (see also |42|). 

A modification of ZM model is proposed by us to adapt 
it for emulation of some aspects of dynamics of microwave 
phonon lasers |2J |6^, 112 (called also phasers 0. HE 



2 



|43() which demonstrate strong instabilities, cooperative 
behaviour, deterministic chaos and other manifestations 
of complexity El The relaxational 
properties of active units (paramagnetic i ons) in phasers 
are of the same type, as in class-B lasers Jsn], and there 
are experimental confirmations |48l . l4f)llH^il53j of common 
properties of phasers and this class of lasers. 

In class-B lasers the inequality T\ ^> Xa e id T 2 
[8l| takes place, where Tfi id is the lifetime of pho- 
tons; Ti and T 2 are respectively the times of longitu- 
dinal and transverse relaxation of active units. The 
same inequality we have for phasers, if one replace pho- 
tons by phonons (see El El 0). From this 
point of view phasers are acoustic analogs of class-B 
lasers. Moreover, acoustic wavelenghts in phasers is 
of the same order as electromagnetic ones in optical 
lasers, because sound velocity in phaser is about 5 or- 
ders less than light velocity in laser (and phasers usu- 
ally operate at frequencies F = ^(P haser ) « l(r 5 F( laser ) 

m, EE E10 E| 00 11 12 ) 

In papers [7a, 1871 l89| class-B lasers were treated 
as a kind of excitable systems, that are widely stud- 
ied in chemical kinetics, especially for the Belousov- 
Zhabotinskii (BZ) reaction and other similar chemical 
phenomena 0, 0, 0, El 0] • Some of differences and 
similarities between class-B lasers and excitable systems 
were examined in [76j: 

"The difference is that excitable systems 
in zero-dimensional case (with spatially ho- 
mogeneous fields) display self-sustained os- 
cillations of the field amplitude, whereas ... 
[class-B systems] display (damped) relaxation 
oscillations only. In the limit 1 711 — » 
the relaxation oscillations decay very slowly; 
thus class-B lasers may be considered quasi- 
excitable systems in this limit. In two spa- 
tial dimensions (unlike the zero-dimensional 
case) class-B lasers display self-sustained os- 
cillations: They behave quite similarly to ex- 
citable s yste ms. The non-stationarity of the 
vortices |83 is one indication of the similar- 
ity with the excitable systems." (0], page 
1143). 

Using pointed similarity and taking into account the 
analogy between phasers and class-B lasers, a phaser dy- 
namics in two-dimension (2D) case may be emulated by 
ZM model El 01 (with some extensions of the last, con- 
cerning relaxational properties of active units). The ZM 
model is a cellular automaton (CA) 1631, i.e. discrete 
mapping, which may be defined as follows. Let an ac- 
tive discretized medium 97l c has the form of rectangular 
2D-lattice. Each cell of the lattice contains one cellular- 
automaton unit (CAU) with coordinates («,j), where 
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(min(i) = 1) A (min(j) = 1). All the CAUs in the 9Jt e are 
identical, and they interact by the same set of rules (CA 
is homogeneous and isotropic in the von Neumann sense 
0). The upgrade of state s[f = S^\i,j) of each CAU 
is carried out synchronously at each step n < N during 
the cellular automaton evolution. The final step N may 
be either predefined or it will be searched during CA evo- 
lution as time (i.e. quantity of discrete steps) for reach- 
ing an attractor of this CA. The conditions of upgrade 
depends both on and on Sj^T , where {«', j'} be- 

longs to certain active neigborhood of CAU at site (i, j). 
In order to formulate boundary conditions for a 9Jt c of 
finite size (i. e. when max(i) = Mx,max(j) = My), a 
set 9Jt v of virtual cells with coordinates (i = 0) V (j = 
0) V (i = M x + 1) V (j = M Y + 1) may be introduced. 

The orig inal ZM model has the single channel of dif- 
fusion 42, 98]. Such one-channel (1C) models are ad- 
equate for chemical reaction-diffusion systems 0]. In 
phaser active system, the multichannel diffusion of spin 
excitations is the typical case, because it proceeds via 
(near)-resonant dipole-dipole (d — d) magnetic interac- 
tions between paramagnetic ions |43|]. For a three- level 
system, which is the simplest phaser system, there are 
3 possible channels of resonant diffusion. In the case 
when d — d interactions are forbidden at one of three res- 
onant frequencies of a three-level system, the asymmetric 
two-channel (2C) diffusion is realized under conditions of 
perfect refractority of the intermediate level — see Ap- 
pendix^ Note that asymmetric diffusion is well known 
in bi ology (see, e.g., the book of D. A. Frank-Kamenetzky 
|l02j ) . Recently a kind of asymme tric d iffusion was pro- 
posed by N. Packard and R. Shaw |103| for a mechanical 
system. 

Another important feature of phaser is extremely low 
level of unavoidable (quantum) noise in its active medium 
because I S pont °c F 3 , where / S pont is intensity of sponta- 
neous emission. Accordingly, Ispont^ ^ s about 15 orders 

lower than because F(P hasor ) /F( lascr ) « 10~ 5 . It 

gives the basis to consider the phaser system as close to 
deterministic, i. e. one need not include stochastic terms 
[85| in CA rules for phaser modeling. 

The paper is organized as follows. In Section[n]we re- 
view deterministic CA generating vortices, giving main 
attention to those of CA, in which dynamically stable 
RSWs are possible. Section IIIII is devoted to detailed 
description of ZM-like TLCA model used in our com- 
puter experiments. The results of computational study of 
RSWs and discussion are presented in Sections llVl andlVl 
Section lvTl contains some concluding remarks. In Appen- 
dices ^ and El the structure of energy levels and relax- 
ation properties of the Ni 2+ : AI2O3 spin-system (used in 
real experiments with phaser) are discussed. Phenomena 
of inversion states collapse and critical slowing down for 
lumped (point-like) phaser models are described in Ap- 
pendix [C] in order to compare them with slowing down 
phenomena in TLCA which is a minimal multiparticle 
model of an active (excitable) medium of phaser. 
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II. ROTATING SPIRAL WAVES IN CELLULAR 
AUTOMATA: A SHORT REVIEW 

Probably the simplest cellular automaton which gen- 
erates rotating (but nonspiral) objects is the Conway's 
Game of Life (CGL) p|. CGL is a two-level (L G £ 
{L , Li}) cellular automaton defined on square grid (97T C 
in CGL is usually infinite) . The state of each cell with co- 
ordinates at step n of evolution is described by very 



simple relationship G = . Here $y ' is binary phase 
counter: £ {0, 1} with such the rules of upgrading: 



(n) 



3> 



(n+i) = 

ij 




where 



(n+1) 



if U% +1) = 2; 
if Ul; +1) = 3; 
otherwise, 



i'ji 



(1) 



(2) 



and the pair of indexes {i' , j'} = {i+a. j+b} describes the 
active neighborhood for the cell For the classical 

CGL this is the Moore neighborhood, i.e. octet of 
the nearest cells: [(|o| < 1) A (|6| < 1)] A (S a0 Sbo / 1), 
where 5 xy is the Kronecker symbol. 

In CGL, rotating (and simultaneously moving) struc- 
tures called gliders are dynamically stable if there are no 
another CGL objects in their vicinity. Gliders are of great 
interest from the mathematical point of view, especially 
because it is possible to construct universal computing 
machine on the basis of gliders and glider guns. On the 
other hand, CGL (at least in its classical two-level form) 
does not describe any known real-world vortices. 

Much more realistic description of vorticity, including 
RSWs, is possible in framework of three-level cellular au- 
tomata (TLCA) or many-level ones. The simplest TLCA 
which generates rotating spiral objects is the single-step 
relaxation (SSR) model, which was extensively investi- 
gated by L. V. Reshodko with co-authors at the begin- 
ning of 1970-th 0,0, H3- Later a similar CA was in- 
dependently proposed and studied by J. Greenberg and 
S. Hastings (GH model) [H- 



A. The single-step relaxations model 

In SSR model, every CAU at each step of evolu- 
tion n occupies one of the three discrete levels Lk £ 
{Li, Lin, ^n}- The order of levels in curly brackets cor- 
responds to the order of CAU's state advancing direction 
(by cycle): Li — > Lm — ► Ljj — ► Lj — ► • • • . Here Li is 
the ground (lowest) level, which is stable in the absence 
of nearest neighbors with K ^ I; Lm is the excited (up- 
per) level, and the Ljj is the refractory (intermediate) 
level. This "arrangement" of levels may be described by 
apparently formal (from the mathematical point of view) 



inequality Lm > Lu > L\. But from the physical point 
of view such the inequality is well defined (see Appendix 
[S]for an example of the concrete physical system). 
The essence of the SSR model is as follows [Tl 121168). 

Let us introduce the phase counters 4>(j <= {~ 1,0, +1}. 
In the SSR model, there is such the correspondence be- 



tween Lk and <p. 



(n). 



(L K = Li) 
[Lk = Lm) 
(L K = L u ) & ( 



= o); 
= +i); 

= -!)■ 



(3) 
(4) 
(5) 



The iterative process (which governs the whole SSR cel- 
lular automaton evolution) describes transitions of every 

individual CAU by mapping of phase counters 



dn). 



a 



(n+l), 
ij \ 



where 



and 



■1, if 4 tt) = +1; 

otherwise, 



(6) 



(7) 



dn) An) \ 



-1, if 



(^=0)A(3* 



(«) 



= +i 



(8) 



0. 



otherwise. 



Here the pair of indexes {i + r,j + s} describes the 
von Neumann neighborhood for the cell (i,j)- This 
neighborhood consists of non-diagonal nearest cells only: 
(|r| = 1) © (| a | = 1) (the symbol ffi means exclusive OR). 

The SSR cellular automaton is one of the simplest 2D 
discrete mappings, which models RSWs in excitable me- 
dia. But the original SSR model has insufficient resources 
to model real physical or chemical excitable systems. In 
particular, transitions in simplest cellular automata of 
the SSR type are "instant" : the lifetimes of excited and 
refractory states are precisely equal to the single step of 
an iterative process. 

Improving of the SSR model may be achieved by 
some different ways: increasing of level numbers ("col- 
ors") and/or extending of the active neighborhood (see 
R. Fisch, J. Gravner, and D. Griffeath [23 for details). 
There are also qualitatively different models which are 
able to generate RSW, as the 2D Oono-Kohmoto CA, 
described in the following Subsection. 



B. The 2D Oono-Kohmoto model 

The well-known ID Oono-Kohmoto model (Y. Oono 
and M. Kohmoto |6^|) may be generalized for the case 
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of 2D active medium (see paper of G. G. Malinetskii and 
M. S. Shakaeva [13 )• Let every CAU has three levels 
{L , Li,Lm}, which corresponds to three discrete values 
of certain global attribute x ( e - g- concentration of a 
chemical agent): € {0,1, M}, where M > 0. The 

evolution of the Oono-Kohmoto 2D (OK2) cellular au- 
tomaton is described by such the discrete mapping: 




1, XW>h 2 ; 

0, if hi < xg° < h 2 ; (9) 



where 

xi n) = (a/AT)X:xg + (l-a)xi n) . (10) 

i'j> 

Here hi,h 2 are lower and upper thresholds for xff 
(e. g. hi = 1/2; h 2 = 3/2), and a characterizes the 
diffusion rate. The pair of indexes {i', j'} = {i +p,j + q} 
describes the active neighborhoods for the cell ftf is 
the coordination number (quantity of cells in the active 
neighborhood). This may be, e. g., the von Neumann 
neighborhood (TV = 4), the Moore neighborhood (TV = 
8) etc. 

G. G. Malinetskii and M. S. Shakaeva revealed RSWs, 
gliders and other interest objects in the OK2 cellular au- 
tomaton [H3 ■ But the OK2 model (as far as the SSR one) 
does not take into consideration relaxation processes. In 
the next Subsection we describe some CA models with 
multi-step relaxation (MSR) mechanisms. 

C. Multi-step relaxation models of excitable media 

Introducing of finite times of relaxation (i.e. delay of 
interlevel transition, which is independent of state of an 
active neighborhood) leads to an important improving 
of cellular-automata models of excitable media. Each 
of non-ground levels in an MSR model has some virtual 
"sublevels" , and relaxation of CAU may be formally de- 
scribed as subsequent jumps between these "sublevels" . 

An extensive investigation of various three-level CA 
with MSR was fulfilled by L. V. Reshodko with co- 
authors in early 1970-th [13, E3, E3l to model excita- 
tions in smooth muscle tissue. The starting point of their 
work was three-level Wiener-Rosenblueth (WR) model 
[9oT |. elaborated in 1946 by N. Wiener and A. Rosen- 
blueth for continuous media. In the original WR model 
|90j . only the refractory level had finite relaxation time. 
The original WR model was modified by A. Rosenblueth 
in 1958 with the purpose to take into consideration 
finite-time relaxation of excited level. The WR model 
was developed in [9(| to model RSWs in inhomogeneous 
media only. But, as it was shown by O. Selfridge [71j and 
I. S. Balakhovskii [13, such the waves are possible in a 
parametrically homogeneous excitable medium too. 



Another successful attempt to use of such the approach 
was undertaken by T. Toffoli and N. Margolus (their re- 
sults, to our knowledge, were firstly published in 1987 
in the book [8(j )• About one year earlier (in 1986) 
V. S. Zykov and A. S. Mikhailov had published paper 
|98| . in which a simple and very clear TLCA was pro- 
posed to model excitable media with arbitrary relaxation 
times (both at excited and refractory levels), arbitrary 
factor of activator accumulation (at ground level) and 
arbitrary active neighborhood. In essence, the model of 
Zykov and Mikhailov j9^ is a generalization of Bogach 
and Reshodko (BR) model [HEl. 

Birth, evolution, interaction and decay of RSWs and 
other spatio-temporal structures in active (excitable) me- 
dia are very sensitive to mechanisms of diffusion. Both 
BR [H[l3| and ZM HHH models are based on the 1C 
diffusion mechanism, which seems to be a good approxi- 
mation for some chemical systems (see, e.g., |83|). Real- 
istic models of several physical excitable system (class-B 
lasers jzl H Hj, phasers HI 0, HI E3 etc.) 
need to take into consideration additional channels of dif- 
fusion, as it was discussed in Section [I] In the next Sec- 
tion we present a detailed description of the MSR model 
of excitable system, based on ZM model and modified 
by us in order to introduce such an additional channel 
in TLCA. This modified ZM cellular automaton with 2C 
diffusion was used in the present paper as the basic mo- 
del for computer studying of RSWs and transient spatio- 
temporal chaos in active (excitable) media. 



III. THREE-LEVEL CELLULAR AUTOMATON 
WITH TWO-CHANNEL DIFFUSION 
(A MODIFIED ZYKOV-MIKHAILOV MODEL) 

A. States of cellular automaton units and branches 
of evolution operator in TLCA 

Let 97t is a rectangular 2D lattice containing Mx x 
My three-level CAUs. The upgrade of state of each 
CAU is (as usually) carried out synchronously during the 
cellular automaton evolution. The excited (upper) level 
L\n has the time of relaxation r e > 1, and the refractory 
(intermediate) level L\\ has the time of relaxation r,. > 1 
(as in the original ZM model [9^ and in contrary to the 
SSR model, where the equality r e = r r = 1 is embedded 
in). Both T e and t t are integer numbers. 

In the TLCA model with 2C diffusion, the first channel 
of diffusion accelerates the transitions L\ — > Lm for a 
given CAU, and the second channel of diffusion (which 
is absent in the ZM model) accelerates the transitions 
Lm — > Lu. The complete description of CAU's states 
in such the TLCA model includes one type of global 
attributes (the phase counters <Py) and two types of 

partial attributes «„■" and z£ for each individual CAU 
in 9Jt . Full description of all CAU's possible states is as 
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follows: correspondences between tpW and Lk take place (by def- 




(11) 



w*n)= T/'(nh : t: (12) 



inition) for all the CAUs in 9Jt e at all steps n of evolution: 

(L K =L 1 )*(<p$ ) =0); (14) 
(L K = Lin) «*■ (0 < < r e ) ; (15) 
(L K = L„) ^ (r e < < r e + r r ) . (16) 



^(L„) = (^). (13) 

In the framework of the TLCA model, the phase coun- 
ters lie in the interval tpffi e [0, r e + r r ]. The following 



These correspondences are of key significance for all 
MSR models (BR model [HEI, ZM model jUi^ etc.): 
there are only three discrete levels, and relaxation of each 
CAU is considered as intralevel transitions (or transitions 
between virtual sublevels). 



The evolution of each individual CAU proceeds by sequential cyclic transitions Lk — > Lk> (where K and K' G 
{I, III, II}), induced by the Kolmogorov evolution operator ft [84] . In the TLCA model, the evolution operator ft has 
three orthogonal branches fti, ftui and fin, which we call ground, excited and refractory branches respectively. The 
choice of the branch at iteration n + 1 is dictated only by the global attribute of the CAU at step n, namely: 



ft 



(n+l) _ 



= ft 



» 



ftm^J, if0<^ 



0, 



< T e 



if T e < < T e 



(17) 



B. Ground branch of the evolution operator 



At step n+l, the branch fti by definition fulfills operations only over those CAUs, which have Lk — L\ at step n. 
These operations are precisely the same, as in ZM model |42,|2g, namely: 




„(«+!) _ A(' 



(19) 



where A^"' is the accumulating term for the it-agent, which is an analog of chemical activator HHi f |" is the 
first-channel diffusion term; h is the threshold for the w-agent (h > 0); g is the accumulation factor for the u- agent 
(g E [0, 1]); C(p, q) is the active neighborhood of the CAU at site and the w-agent arrives to CAUs with Lk = L\ 

only from CAUs with Lk = im in C (p, q) : 



J, 



(n) 

i+p,j+q 



1, if (0<^; j+9 <r e ); 
.0, if (vffi, li+9 >re)v(^ ii+ , = 0). 



(20) 



The definition of the diffusion term -D l -™' ) in Eqn. i|19|) is very flexible. Apart of well-known neighborhoods of the 
Moore (C(p, q) = Cm(p, <?)) and the von Neumann (C(p, q) — Cn(p, q)) types: 



Cu{p,q) = 



1, if 



( N < 1) A ( |g| < 1) A (<5 p0 <V ? 1) 



0, otherwise, 



(21) 
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CM J 1 ' *[(IH-i)«(M-i)]; (22) 

[0, otherwise, 

one can easy define another neighborhoods, e. g. of the box type Cb(p, q) or of the diamond type Cd(p, q) |I3 (which 
are the straightforward generalization of the Moore Cm{p, q) and the von Neumann Cn{p, q) neighborhoods), etc. 

On the other hand, the ZM definition 0,113] of the weight factors jj-™*, - + if^Hjl may be extended out of the binary 
set {0, 1} to model various distance-dependent phenomena within, e. g., C 'd(p, q), C 'b(p, q) ■ 

In this work, however, we restricted ourselves by the Moore neighborhood and by weight factors of the form 
(|20p. Another types of neighborhoods, weight factors and some other modifications of the model will be studied in 
subsequent papers. 



C. Excited branch of the evolution operator 



At the same step n+1, the branch flm fulfills operations only over those CAUs, which have Lk = Lm at step n: 



in 



(*o 

<Pi3 




if 



(0 < ^ < r e ) A (z^ < f) V = r e ) 



y ("+l) 



if (0 < < T e) A > /) 



(23) 



4? +1) = 4 B) - E c fe 9)qS j+9 . (24) 
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where is the second-channel diffusion term; / is the threshold for the z-agent (/ > 0), and we assume that 

z-agent arrives to excited CAU (Lk = Lm) only from those CAUs, which have Lk = L\ in C(p, q): 



Q ^-\o, if^L^o. (25) 

One can see from (|24|l that the z-agent does not accumulate during successive iterations. In other words, the branch 
Sim at step n+1 produces "memoryless" values of partial attributes zj™ for CAUs having Lk = Lm at step n (in 
contrary to the u-agent for CAUs having Lk = L\ at step n). The z-agent may accelerate transitions from excited 
CAUs to refractory ones. This is important difference between 2C model and the original 1C model of ZM [42l l98| : 
and 2C model qualitatively describes the real active media used in quantum acoustics (see Appendices 1X1 and iBl) . 



D. Refractory branch of the evolution operator 



The branch fin does not produce/change any partial attributes at all (because the intermediate level L\\ is in the 
state of refractority). It fulfils only such the operations over CAUs, which have Lk = L\\ at step n: 

V J (0, if (flj' = T e +T r . 

Generally speaking, there are many examples of active media with weak refractority (when the unit at the inter- 
mediate level L\\ is not absolutely isolated from its neighbor hooding units). But in this work we restrict ourselves by 
the case of perfect refractority l|26|l . which is valid, e. g., for the phaser systems of the Ni 2+ : AI2O3 type (as it was 
pointed out in Chapter [I] and is explained in Appendix |A"}| . 



E. Reduction of the TLCA model to the ZM automaton 



The second- channel diffusion gives the contribution to the TLCA dynamics if / < Af, i.e. if / < 4 for C = Cjv, 
/ < 8 for C = Cm and so on. 
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At ((C = Cat) A (/ > 4)) our TLCA model is of ZM-like (i.e. 1C) type, and at ((C = C M ) A (/ > 8)) it becomes 



equivalent to the original ZM model 42, 98], which (with slight rearrangement of cases) is as follows: 



(n+l) 

<Pii = < 



(n) 

1, 
0. 



if < tptf < r e 



•r ( (n) 

if w = 

» 



0)A(u^>>h)- 1 

(n+l) 



W = 0) A (tig* 1 ' < fc) VG#>=T. + T P ), 



(27) 



where m-™^ is defined as in ljT§|) ~ (|2~T]l . 



F. Comparison with other models 



In our computer experiments we used such the condi- 

tions: G I ! 1 }) A e i ^})- Under 

these conditions and in particular (but very important) 
case of g € {0; 1}, the ZM model HlHi becomes of fully 
integer kind. If all these conditions are supplemented by 
Q[+ P , J+q G {0; 1}, the TLCA model Q)-© becomes 
fully integer too. Finally, at C = Cm, 9 = 0, h = 1, / > 8 
the TLCA model I|11[) - (I26|I is equivalent to the simplest 
variant of the ZM modelj42l l9S| , which corresponds to a 
discrete form of the original WR model [9(| . 

On the other hand, the TLCA model l(TT )l -l(25 )l (and, 
usually, the original ZM model |42l 198^ may be consid- 
ered as the generalization of the SSR model. This gener- 
alization is different from those of R. Fisch e. a. ,22] in 
many aspects; first of them consists in using a concept of 
multilevel transitions in [2^] instead of a concept of re- 
laxation times, used in works |42ll98| and in the present 
work, formulae (|I1|I to (|26l) . But in the limiting case of 
"instant" (i.e. one-step) transitions both TLCA model 
and the model of R. Fisch e. a. [22| reduce to the same 
original SSR model. For the TLCA model this is the 
case of C = C^,f > A,g = 0, h = l,r e = r r = 1. 
Under these conditions the system of equations Qllfl- 

©-© if onc uses tne ansatz (^>-™ € 

b\f 6 {2,0, 1}) in the SSR model. Such 
the ansatz is correct because the SSR model in the form 
of Js3 - © is almost insensitive to the phase counter value 
for Lif. the "—1" in the equations (JSJ and Q may be 
changed to any X £ {0, +1 }. 

At (r e > 1) V(tv > 1) our TLCA model is of MSR type, 
and relaxation of each individual CAU is described in the 
same manner as in various Reshodko's MSR models of 
excitable media but diffusion mechanism in TLCA 
model is algorithmizcd using more versatile way. 

There is no direct correspondence between TLCA 
(ZM) model and the OK2 model introduced in de- 
spite of definite similarities in spatio-temporal dynamics 
of these CA. Partial correspondence is between TLCA 
(ZM) model and one of variants of Toffoli-Margolus mo- 
del, namely with its "monotonic-threshold" variant (see 
Chap. 9 in H3). The TLCA model reduces to that vari- 
ant of Toffoli-Margolus model if C = Cm,/ > 8,5 = 
0,he {2,3,4},r e = l, Tr = 2. 



(|2*S|l reduces to 

{-1,0,1})- 



G. Geometry of active media, boundary conditions 
and transients in CA 

Pattern evolution is very sensitive to geometry of active 
medium and boundary conditions even in the framework 
of operation of the same CA. There are three types of 
geometry and boundary condition (GBC) combinations 
most commonly used in computer experiments with CA 
(see Table I). Let us consider them in details. 

• The GBC-1 type is defined as follows: 



where 




if i v = 0; 

if i v = M x + 1; 

otherwise, 



if j v = 0; 

if j v = M Y + 1; 

otherwise. 



(28) 



(29) 



(30) 



Here (i v ,j v ) £ 9Jt v , and the coordinates (ig,je) belong to 
the border set QTlJT C Wl c of excitable area, i. e. (i e ,je) — 
(C j c r ) if (i = 1) V (j = 1) V (i = Mx) V (j = M Y ). ' 

• The GBC-2 type can be easily defined by using an 
extended interval for the phase counter (p. Let every 
cell (i v ,j v ) in 97t v contains one unexcitable CAU. All 
these unexcitable CAUs have frozen (p^ n '(i v ,j v ) a t all 
n E [0,N] without reference to states of («e,if), and 



Jv) = T e 



(31) 



where, a > 0, say a = 1 without loss of generality. These 
unexcitable CAUs in 97l v "absorb" excitations at the bor- 
der of 33T e in contrary to the case of GBC-1, where excita- 
tions at the border 3Jl are reinjected in active medium. 
In many cases this difference leads to qualitatively differ- 
ent behaviour of the whole CA. Note that by this way one 
may also define CA with inhomogeneous active medium 
(where some of CAUs are unexcitable not only in SDt v , 
but in 5DT too), introducing an additional orthogonal 
branch fig into the evolution operator Q. This branch 



is activated at ip. 



T e + T r + 1 only and is simply the 
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TABLE I: Geometry of active media and boundary conditions 



Type 


Geometry 


Boundary conditions 


GBC-1 


Bounded, Toroidal 


Cyclic 


GBC-2 


Bounded, Flat 


Zero-flow 


GBC-3 


Unbounded, Flat 


(free space) 



between the TLCA model and the real phaser systems 
|1 IH |H © |H |H |H HI HI, but in any case au- 
tonomous phaser has only these two mechanisms of in- 
teraction with the outer world. 



H. Initial conditions 



identity operator £1 = /, i.e.: y>-™ +1 ' = fioV« = <Pij ■ 
And even more complex behaviour of such "impurity" 
CAUs may be defined using analogous approach (i. e. by 
introducing additional orthogonal branches in an evolu- 
tion operator): there may be pacemakers |42j or another 
special units embedded in active medium and described 
by phase counter in extended areas (ip > r e +T r )\/(ip < 0). 

• The GBC-3 type is, strictly speaking, a case of bor- 
derless system. But the starting pattern has, of course, 
bounded quantity of CAUs with Lk Lj, located in 
bounded part of active medium. Hence, during the whole 
evolution (1 < n < N) the front of growing excited area 
will meet the unexcited (but excitable!) CAUs only. For 
the case GBC-3, there is neither absorption of excitations 
at boundaries (as for GBC-2), nor feedback by reinjec- 
tion of excitations in the active system (as for GBC-1). 
From this point of view CAs of GBC-3 type are "sim- 
pler", than CAs of GBC-1 and GBC-2 types. On the 
other hand, CAs of GBC-3 type are potentially infinite 
discrete systems where true aperiodic (irregular, chaotic) 
motions are possible, in contrast to finite discrete sys- 
tems, possessing, of course, only periodic trajectories in 
certain phase space after ending of transient stage .41j . 

On the other hand, there is also very important dif- 
ference between GBC-2 and other two types of GBC. 
The system of GBC-2 type is the single from the three 
ones under consideration which interacts with the ex- 
ternal world dynamically (besides of relaxation). Sure 
enough, GBC-1 and GBC-3 are connected to this outer 
world only by relaxation channels (the r e and r r are the 
measures of this connection). In contrary, a CA of the 
GBC-2 type interacts with surroundings through the real 
boundary, which is in fact absent in toroidal finite-size 
active medium of GBC-1 type and it is absent by def- 
inition for CA of GBC-3 type. Despite of elementary 
mode of such interaction (boundary simply "absorbs" 
the outside-directed flow of excitations from a 
CA of the GBC-2 type may demonstrate very special 
behaviour. The main attention in this work is devoted 
just to TLCA of GBC-2 type as the most realistic mo- 
del of phaser active system, where both the mechanisms 
of interaction (dynamical and relaxational) of an active 
medium with the outer world are essential. Dissipation in 
a phaser active medium (highly perfect single crystal at 
liquid helium temperatures) is caused by two main mech- 
anisms: (a) dynamical, by coherent microwave phonon 
and photon emission directly through crystal boundaries 
and (b) relaxational, by thermal phonon emission. There 
are, of course, many more or less important differences 



Patterns P'™- 1 in CA are usually described in terms of 
levels (or "colors") of CAUs, and initial conditions P^ ' 
are formulated simply as matrix of CAUs levels Lk- But 
states of CAUs in such the automata as ZM or TLCA 
are fully defined not only by levels Lk itself. There are 
global attributes which must be predefined before TLCA 
evolution is started. Some of parial attributes must be 
predefined too. These points are essential for reproducing 
of the results of computer experiments. 

In our work such the initial conditions for the global 
attributes ipij are used: 

<4° ) = 1 > XLW=L m ; (32) 

^r. + l, ULf=L lh 

Initial conditions for mj must be defined for ground-state 
CAU's (Lj( = L\) only. In this work, we suppose 

(4 0) =0) IFF (<p<g>=0), (33) 

where IFF means "if and only if" . Initial conditions for 
Zij are undefined for all Lg because z-agent is not de- 
fined at n — 0, see Eqn. (|12|) . As the result, starting 
pattern P<°) may be defined (both for the ZM and TLCA 

models) as the matrix ||<£-° || (where <p^ £ {0, l,r e + l}) 
with additional condition given by Eq. I|33[) . 

IV. RESULTS OF COMPUTER EXPERIMENTS 

Cellular automata are inherently based on irreducible 
algorithms. Generally speaking, there are no predictive 
procedures for such the systems. So, the best way to 
investigate cellular automaton (TLCA in particular) is 
to run it, because, as S. Wolfram pointed out, "their own 
evolution is effectively the most efficient procedure for 
determining their future" (see j92, page 737). 

Here we present results of our computer experiments 
both with original 1C automaton of ZM and with 2C 
automaton described in Section IIIII The main tool in 
these experiments was the program "Three-Level Cellu- 
lar Automaton" (TLCA © 2004 S. D. Makovetskiy 
which is included now in the software package "Three- 
Level Laser Model" 99]. Original ZM-algorithm HHHl, 
was implemented in a previous software package "Gen- 
eralized Wiener-Rosenblueth Model" (GWR © 2001- 
2003 S. D. Makovetskiy [H; see also [34^. 

Computer experiments were fulfilled for wide ranges of 
control parameters (CP): (1 < r e < 50) A (1 < r r < 50) 
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(for all the cases); < g < 1 (for non-integer versions 
of TLCA and GWR); g £ {0; 1} (for integer versions 
of TLCA and GWR); 2 < h < 100 (for all the cases); 
1 < / < 9 (the case / = 9 corresponds to original ZM 
model — see Subsection IIII E|) . 

A discrete system with finite set of levels may have only 
two types of dynamically stable states at bounded lattice. 
The first of them is stationary state, and the second is 
periodic one. They may be called attractor s by an alog y 
with lumped dynamical systems (see e. g. |4(1 l93t 19^^. 
For our cellular automaton, the first of such the attrac- 
tors is spatially-uniform and time-independent state with 
<4 l) = 0, where (i € [1,M X ]) A (j e [1,M Y ]), n > n c . 
In other words, the only stationary state of TLCA is the 
state of full collapse of excitations at some step nc (by 
definition of excitable system) . The second type of TLCA 
attractors includes many various periodically repeated 
states (RSWs is a typical but not the single case) . In this 
case V f = 4" +T) where {% £ [1,M X ]) A (j £ [1,M Y ]), 
n > np] np is the first step of motion at a periodic at- 
tractor; T is the integer-number period of this motion, 
T > 1. 

If starting patterns (where ID is the pattern 

identifier) are generated with random spatial distribu- 
tion by levels, then the time intervals nc or np may 
be considered as times of full ordering in the system. 
Such irregular transients are of special interest from 
the point of view of nonlinear dynamics of distributed 
systems B El US H ilj because they may be bot- 
tlenecked by very slow, intermittent morphogenesis of 
spatial-temporal structures. In the next Subsection we 
study this collective relaxation both for cases of collapse 
and periodic final states of TLCA evolution. 

A. Bottlenecked collective relaxation in TLCA 

Evolution of TLCA usually has several stages with very 
different characteristic times and qualitatively different 
spatio-temporal dynamics. An example of relatively sim- 
ple evolution is shown at Figure 1 for starting pattern 
p (o) _ p (o)_2 Here j = 8j but at / = 9 this starting pat- 
tern reaches precisely the same attractor. This is because 
our TLCA model at / = 8 is very close to ZM model. 
But even at / = 8 such coinciding of attractors is rela- 
tively rare (less than 10 percent for the CP set pointed 
at caption to Figure 1 and starting patterns of the type 
shown at Figure Al). In other words TLCA model even 
at / = 8 is not equal to ZM model. Much more typical 
case is gradual divergence of spatio-temporal structures 
for/ = 8 and / = 9 during evolution (see Figure 2). And 
for / < 8 we have not observed any attractors coinciding 
with ones for / = 9. 



Evolution of TLCA with / < 8 is more complicated, 
and transient time may rich giant values (millions of steps 
for the same Mx = My = 100 as for patterns at Fig- 
ure 1 and Figure 2) due to bottlenecked RSWs mor- 
phogenesis. Under certain conditions, the morphogen- 
esis of RSWs proceeds by multiple irregular changings 
of effective topolgical charge 3 Qt, including reversing of 
sgn(QT)- For example, at / = 4 the system with starting 
pattern P(°) = P^ (see the P^ at Figure Al) reaches 
its periodic attractor only at n ~ 4 ■ 10 6 . Seemingly per- 
fect RSW with Q T = +2 is forming at n « 3 • 10 5 , but 
it is obviously not an attractor. At n ~ 4.5 ■ 10 5 the 
core of RSW becomes complicated and gradually evolves 
to state with Qt — +3 (Figure 3, step n = 500525). 
During sequence of several metamorphoses, the system 
returns to Qt = +2 (Figure 3, step n = 1.35 ■ 10 6 ). And 
at n = 1.63 • 10 6 one can see RSW with opposite sign 
(Qt — —2). The further evolution continues this un- 
predictable scenario, ending at attractor with Qt = —3 
(steps n = 3.97- 10 6 to n = 4- 10 6 at Figure 3). Besides of 
quantity and sign of Qt, the frequency of spatial waves 
irrregularly changes too — compare, e. g., patterns at 
n = 2 • 10 5 and n = 1.09 • 10 6 (or at n = 1.94 • 10 6 and 
n = 2.59 • 10 6 etc.). 

It is interesting to compare such scenario of evolution 
of TLCA with a scenario of crystal growth due to spiral 
(screw) defects increasing. The results are just opposite: 
crystal growth is usually accelerated by many orders due 
to screw defects, but TLCA relaxation (transient to an 
attractor) is highly bottlenecked by RSWs. 

There is a huge amount of periodic attractors for such 
TLCA, but the steady-state attractor is the single one. 
It is fully collapsed state of an excitable system, as it 
was pointed out earlier. On the other hand, there is 
very many starting patterns that evolve to this collapsed 
state (our TLCA is an irreversible system). An example 
of very slow transition to collapse is shown at Figure 4. 

Bottlenecked collective relaxation shown at Figure 3 
and Figure 4 is, of course, aperiodic. But it has some 
quasiperiodic features because observed at Figure 3 and 
Figure 4 spirals (or spiral "domains" ) possess more or less 
regular structure. Much more chaotic spatio-temporal 
behaviour one can observe varying CP of the TLCA. Such 
a complex, slow and unpredictable evolution of a simple 
TLCA system (with / = 9) is shown at Figure 5. This is 
a typical example of transient spatio-temporal chaos in 
collective relaxation of excitations, which we observed in 
TLCA. 

Such highly bottlenecked and irregular transient dy- 
namics (metaphorically, "turbulence" ) was observed nu- 
merically by J. P. Crutchfield and K. Kaneko 18] in an 
ID model which is a kind of coupled map lattice (CML) 
model 32]. In contrary to CA and other discrete map- 



2 Typical set of starting patterns with Mx = My = 100, used in 
our numerical experiments, see Figure Al. 



3 Effective topolgical charge is defined in the near vicinity of the 
vortex core. 
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pings, CML is not fully discrete system because of contin- 
uous spectrum of states of its elementary units. In other 
words, CML has not any well resolved levels. But both 
TLCA and CML demonstrate lethargic transients which 
are not consequences of usual critical slowing down (i. c. 
singularity of transient time at several combination of 
CP [23 ) ■ Critical slowing down takes place near critical 
points only, and fine tuning is needed to reach the (very 
narrow) band of CP where this type of slowing-down is 
observable. 

In phaser amplifier the observed phenomenon of crit- 
ical slowing down |m| may be described by simple 
lumped model of inversion states, see Appendix [CJ But 
in autonomous phaser generator we experimentally ob- 
served [44j phenomena of spin-phonon interaction with 
very long, lethargic aperiodic transients (phonon "tur- 
bulence" ) which have qualitatively another nature than 
usual critical slowing down. 

There are some qualitatively different mechanisms of 
slowing down, for which such the tuning is not necessary. 
A well known mechanism of this type is the effect of 
self-organized criticality ;9|, which takes place when a 
dissipative system holds themself in critical state and no 
external tuning of CP is needed (models of sandpiles, 
earthquakes, forest fire etc.). 

Some interesting scenarios of bottlcnccked relaxation 
we re p roposed recently by K. Kaneko and his co-authors 
particular, self-organized bottleneck was 
revealed by H. Morita and K. Kaneko [6(j in transient 
processes for simple excited Hamiltonian system. Here 
"the critical state is spontaneously formed without con- 
tinuous driving, once a part of a system is highly excited" 
60]. And in work of A. Awazu and K. Kaneko Q it 
was found that forming and evolution of transient dis- 
sipative structures lead to non-critical slowing down in 
closed chemical system (a modified Brusselator model). 

Super-slow evolution of seemingly "frozen" spiral 2D 
structures was revealed and investigated by another 
group (C. Brito, I. S. Aranson, and H. Chate) pjj) m 
2D system modeled by complex Ginzburg-Landau equa- 
tion. During last years some new publications on slow 
evolving continuous systems has been appeared, partic- 
ularly in the "Condensed Matter" and "Nonlinear Sci- 
ences" arXives. 

Collective relaxation phenomena observed by us in 
TLCA are of similar kind despite of fully dicrete nature 
of our model. Bottlcnccked transition to final state of 
TLCA caused by RSWs may qualitatively explain why 
relaxation times of level-populations in dissipative three- 
level (or multilevel) paramagnetic systems differ by sev- 
eral orders from those predicted by one-particle spin- 
lattice relaxation theory 0, [2(j (this old problem has 
not satisfactory solution up to now). On the other hand, 
freezing, quasi-freezing and slow evolution are now a sub- 
ject of investigation of three-level systems in many alter- 
native directions of modern nonlinear science, see e . g . 
work of F. Vazquez, P. L. Krapivsky, and S. Redner [8q 
etc. 



Hypersensitivity of the TLCA to changing of initial 
conditions (which will be discussed in the next Subsec- 
tions of this work) , instability of vortices and extremely 
slow, lethargic moving to a regular attractor are typical 
features of transient deterministic chaotic behaviour of 
TLCA. Birth, evolution, interaction and decay of RSWs 
in TLCA (including original ZM model), as it is obvi- 
ous from Figures 1 to 5, is much more complex and un- 
predictable than it was expected at early days of inves- 
tigation of such axiomatic CA models of excitable sys- 
tems. In particular, collisions of RSWs with absorbing 
borders (in the GBC-2 case) violate conservation of Qt 
(see Figure 3 and Figure 5). In the next Subsection we de- 
scribe some our numerical results concerning the problem 
of competition of left- and right-handed spatio-temporal 
structures, which is of interest not only from physical 
point of view |l5j . 



B. Competition of left- and right-handed RSWs 

A bounded solitary domain of excited CAUs having 
appropriate relaxation times and placed far from grid 
boundaries (or at unbounded grid) may evolve to RSWs 
if and only if there is an adjoined (but not a surround- 
ing) domain of refractory CAUs [lfl |H |H [FJ HjJ . In 
this case RSWs apear by pair with opposite sgn(QT) and 

integral Q T ntesT ^ = J2 m Qt i s obviously conserved (by 
infinity in time if grid is unbounded). If such solitary 
excited-refractory area is placed in a starting pattern 
near the grid boundary and the GBC-2 conditions take 



place, the single RSW appears and evolves and Q T nteBT ^ 
is, of course, not conserved in this case. These simplest 
and well known examples illustrate possibilities of coex- 
istence and competition of one or two RSWs in excitable 
media, and possible scenarios of the RSW(s) evolution 
may be easily forecasted. 

An evolution of complex patterns with multiple, irreg- 
ularly appearing, chaotically-like drifting and colliding 
RSWs (or, possibly, another spatio-temporal structures 
posessing handedness) is in essence umpredictable with- 
out direct computing of the whole transient stage. In 
particular, some intuitive predictions about of Q^ ntesr ^ 
at attractor may be fully errorneous, as it will be shown 
in this Subsection. 

In the previous Subsection, it was already pointed out 
that, for some CP sets, an effective topological charge Qt 
(more precisely, Q T nteBT ^) irregularly changes, including 
reversing of its sign during a system evolution. Extremely 
high complexity of almost all pre-attractor evolution (see 
e. g. Figure 5) does not permit to prognose even the 
near future of such TLCA. But sometimes it seems that 
long-term prediction is yet possible for another CP sets. 
In particular, Figure 6 demonstrate such "possibility" at 
least at pre-finish stage of transient process in TLCA 
with / = 3 (other CP are the same as for Figures 1 
to 4) and with Mx = My — 300. Indeed, beginning 
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from n > 10 5 the sign of the integral effective topological 
charge becomes equal to sgn ^Q^ ntcgr ^ = — 1 (Figure 6, 

n = 1.6 • 10 5 ) and stays unchanged up to winning of the 
single RSW in the spiral competitions. In other word, a 
"perfect" left-handed word is formed here. 

Nevertheless, this pretty picture of "predictable" self- 
organization is fully destroyed at n « 3.96 • 10 5 due to 
collision of the winner's core with the grid boundary and 
subsequent full collapse of excitations at np = 396982. 
This last stage is not shown at Figure 6 because it is of 
the same type as collapse at Figures 4 and 5 - it is simply 
a gradual absorbing of non-spiral waves of excitation by 
boundaries of our GBC-2 system. 

But much more unexpected, fully counter-intuitive 
TLCA dynamics may be illustrated by two snapshots (see 
Figures 7 and 8) of evolution of a system with the same 
CP set as for Figure 6, but with M x = My = 900. At 
Figure 7, one can see almost purely right-handed pat- 
tern at step n = 1.48 • 10 6 of TLCA evolution. The 
only deviation from ideal right-handedness at Figures 7 
is the rudimentary counter-clockwise "tail" of the clock- 
wise spiral at top-center part of Figure 7 (note that this 
highly asymmetric spiral- antispiral pair is inside of the 
closed wavefront). 

One may expect that the winner RSW in this system 
will be of the same handedness (by analogy with the Fig- 
ure 6, where left-handed multi-spiral pattern was evolved 
to left-handed solitary RSW) . But the subsequent evolu- 
tion of the system is very surprising, see Figure 8. The 
winner RSW is shown here at step n = 2 • 10 6 , it rotates 
counter-clockwise {Qt = — 1), i- e. the system becomes 
left-handed despite of their long right-handed previous 
life. During further evolution [n > 2 ■ 10 6 , not shown 
here), the winner's core is moving irregularly across the 
grid up to collision with the grid boundary, and full col- 
lapse of excitations takes place at nc — 2260964 with 
final stage by scenario shown at Figures 4 and 5. 

We will emphasise that the single left-handed RSW 
(Figure 8) is the result of evolution of the right-handed 
pattern of RSWs (Figure 7). A series of intermediate 
snapshots (not shown here due to their big size) at the 
time interval from n = 1.48 • 10 6 to n = 2 • 10 6 demon- 
strate a scenario of such reversing of the system handed- 
ness. An the cause of this reversing is not rudimentary 
counter-clockwise "tail" of the clockwise spiral at top- 
center part of Figure 7, because it quickly dies during 
"spiral war" . Moreover, at n = 1.88 • 10 6 the single fully 
developed right-handed RSW remains at the whole grid, 
in the right-bottom corner of the grid. The second com- 
peting right-handed RSW is in a dangerous proximity 
from left border of the grid at the same n = 1.88 • 10 6 . 
And a small dislocation in spiral wavefront of the first 
RSW appears as the nucleus of future left-handed RSW. 

At n = 1.89 • 10 6 , the two competing right-handed 
RSWs are still alive (both of them are near corners of 
grid), but the mentioned dislocation has already tran- 
sormed to fully developed left-handed RSW. During the 



subsequent evolution, both right-handed RSWs are died, 
the remnant of wavefronts of right-handed structure is 
crowded out by giant lacuna, and the left-handed RSW 
becomes the winner in the world with almost perfect 
right-handed history. Finally, this solitary left-handed 
RSW must die too at nc ~ 2.26 • 10 6 , as it was pointed 
already 

This intriguing story confirms high level of unpre- 
dictability of the system under consideration, and this 
subject will be discussed in details below (see Section 
IVj) . And now we describe another unusual form of mo- 
tions in our TLCA observed at different values / = 5 of 
the threshold for the 



C. Chimera states - coexistence of regular and 
chaotic domains in TLCA 

At the same CP T e ,r r ,h,g as in Figures 1 to 4 and 
Figures 6 to 8, but with / = 5 we revealed unusual, long- 
living stages of CA evolution, during which coexistence of 
periodic and aperiodic spatio-temporal structures takes 
place in TLCA. A useful tool for investigation of such 
the phenomenon is generalized Poincare cross-section of 
pattern sequence. 

Ascending and descending generalized Poincare cross- 
sections, used in this work, are defined as follows. Let 
B(n) is an 1-bit (e. g. black = TRUE, white = FALSE) 
pattern and An is a selection interval (An > 1). For 
each binary cell By 6 B, wc fulfill the sequence of k 
logical additions in the ascending time (starting from a 
step n\ < N — kAn): 

B.[+ ) (n 1 ,An,fc)=[[[B y (n 1 )V 

VB^(ni + An)]V (34) 
VB^(ni +2An)] V ... 
WBij(ni + kAn) ], 

and a sequence of m logical additions in the descending 
time (starting from a step n 2 > mAn): 

B[-\n 2 ,An,m) =[ [ [ By(n 2 )V 

VBy(n 2 - An)]V ^ 

V By(n 2 - 2An) ] V ... 

V By(n 2 — mAn) ]. 

Ascending and descending generalized Poincare cross- 
sections (V^ and V^) are the patterns, which consist 
of B\j(ni, An, k) and B,[~' ) (n 2 , An, m) respectively: 

pW(m,An,fe) = [|B^ +) (m,An,fc)|| (36) 



V^(n 2 ,An,m) = ||B t ^(n 2 , An,m)|| (37) 



12 



Figure 9 demonstrates phenomenon of coexistence of 
long-living periodic and aperiodic spatio-temporal struc- 
tures in TLCA by using of ascending and descending gen- 
eralized Poincare cross-sections (left and right columns of 
Figure 9 respectively; the central column at Figure 9 is 
usual sequence of TLCA patterns showed for clarity) . We 
use term "transient chimera states" for such structures 
by analogy to term introduced by D. M. Abrams and 
S. H. Strogatz [j| for analogous non- transient phenom- 
ena revealed earlier by Y. Kuramoto and D. Battogtokh 
[38| (see below in this Subsection). 

The series of ascending generalized Poincare cross- 
sections V^ + \n,x, An, k p ) shows that spatial domain of 
aperiodic motions is bounded at the whole interval of 
selection. In other words, periodic domain is dynami- 
cally stable (beginning at list from 1.5 • 10 5 ) for more 
than half of the grid area. And the series of descend- 
ing generalized Poincare cross-sections 7- , ' _ * ) (ri2, Art, m q ) 
demonstrates growing of the regular domain, i. e. how the 
system reaches the attractor keeping previously occupied 
territory of local periodicity. This evolving coexistence 
of two qualitatively different (periodically and aperiod- 
ically oscillating) domains in high-dimensional discrete 
mapping we call transient chimera states. 

For r e = r r = h = 50, .9 = 1 and / =/= 5 reaching of an 
attractor proceeds typically by metamorphoses of aperi- 
odic waves to periodic ones over the whole grid, without 
coexistence of periodic and aperiodic domains. But for 
some alternative values of r ei T ri h, g and for / =/= 5 we ob- 
served chimera states again. Evolution of some of these 
states are different comparatively to ones shown at Fig- 
ure 9. For example, small but strictly regular spatial 
islands in chaotic sea may appear; they live by long time 
and then dissappear during that or those stages of TLCA 
evolution; and attractor may be reached after disappear- 
ing of all of these islands. In other words, emergence of 
regular domain(s) of CAUs (shown at Figure 9) is not the 
single scenario of an attractor reaching when spatially- 
periodic and spatially-chaotic states are coexisting. 

Similar phenomena of coexistence of incongruous or 
even antagonistic spatio-temporal structures were dis- 
covered during last years in computer experiments with 
ID arrays of identical coupled oscillators 0, |3^ and 2D 
grids emulating blackouts in electric power systems [l7|. 
The most essential difference between systems, which was 
used in [si UtI Issl ] and in the present work, consists in the 
nature of their elementary units. In works [3J, l38j and in 
overview |17| one deals with units, which by definition 
oscillate even without any coupling with their neighbor- 
hoods. In the TLCA model (as well as in all models of 
excitable media) each isolated elementary unit has the 
single stable state, namely steady ground state. 

In other words, the TLCA or another system of this 
kind may be treated as oscillator if (and only if) there is 
more or less strong interaction between elementary units. 
And namely the whole TLCA is an oscillator in this case. 
Electric power generators or another monochromatic os- 
cillators, of course, may work autonomously — with their 



own frequencies and phases. But a grid with noninter- 
acting excitable CAUs at n > r e + r r is simply a collec- 
tion of anabiotic cells. So the concepts of phase synchro- 
nization/desynchronization are unapplicable to individ- 
ual CAUs by definition. But synchronization may take 
place between domains of CAUs (and this is not a triv- 
ial "synchronization" of e.g. pacemaker-like domains in 
excitable media with stable domain walls) . 

To conclude this Subsection, we mention that coexis- 
tence of static regular and disordered structures is well 
known in CGL (complex "still life" configurations of 
CAUs). E. g., some unbounded (GBC-3 conditions) two- 
level (Lg G {Lq.Li}) CGL-like automata demonstrate 
unlimited growing of domain(s) with frozen CAU states 
(Turing structures), were static spatially-chaotic "cloud" 
is irregularly grided by static spatially-periodic "rays" . 
The simplest CA of this type is the INKSPOT automa- 
ton which may be defined similarly to CGL automa- 
ton (see Eqns. (QJ and @), using binary phase counter 
G {0, 1} and the following elementary upgrading 

rules: 



(n+l) 



1, 



(») 



if U^ +1) + 3; 

if ui; +1) = 3, 



(38) 



where U^ +1 ^ is defined by Eqn. J3J) and the Moore neigh- 
borhood is used. 

Only Lq — > L± transitions are permitted in the 
INKSPOT automaton (|38|l . but this restriction is not 
critical for the phenomenon of coexisting of regular and 
chaotic static domains. A similar behavior demonstrate, 
e. g.. the RtC ("Rays through Clouds" ) two-level automa- 
ton 46] defined by rules 



(n+l) 




if (^ +1) = 2 )v(l^ +1 >G[4;8]) ; 



ift# +1) =3; 
otherwise, 



or, in more convenient form: 
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(n+l) 
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(n+l) 

€ {0; 1}, 



{0;1;3}; 
3; 



(39) 



(40) 



It is clear from the equation (|40|l that both transitions 
Lq <=^ L\ are permitted for the RtC automaton, but 
emergent Turing structures of chimera type is still the 
main form of evolution of this C A (except of special start- 
ing patterns giving pure rays) |4r| . 



D. Spatial symmetry breaking and restoring 

Computer experiments with spatially random input 
leave out an another important case of initialization of a 
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system by symmetric or almost symmetric starting pat- 
tern. In this Subsection wc investigate TLCA evolution 
with such (almost) symmetric initial conditions. 

For a perfectly symmetric input, the TLCA evolu- 
tion will proceed by conserved-symmetry scenario. It 
directly follows from the TLCA definition as homoge- 
neous and isotropic (in the von Neumann sense [62|) Kol- 
mogorov machine with deterministic rules. But how will 
the TLCA evolve if a perfect symmetry of input pattern 
is more or less distorted? Will symmetry be broken out or 
it may be restored? Or both these scenarios are possible 
due to multistability of TLCA? One cannot find answers 
to these questions a priori, using only general guidelines 
as Curie's principle fl4l |. 

To find such the answers empirically, let us study 
TLCA evolution with strictly determined asymmetry of 
staring patterns (Figure 10). 

Left-column starting pattern P$l at Figure 10 is not 
fully symmetric (ten black cells at the top of the excited 
left part of P^jl are changed to ten white cells at the top 
of the excited right part Pml)- But the C v symmetry is 
fully restored already at n = 48 (not shown at Figure 10), 
and at n = 10129 the system is already at periodic at- 
tractor ("Smile of Cheshire cat" 4 — see Figure 10, left 
column, row d). This attractor has period T = 3592 and 
consists of two mirror-symmetric, counter-rotating RSW 
with \Q T \ = 3. 

The right-column starting pattern Pmr at Figure 10 
differs from the left-column starting pattern P^-^ only 
by the single CAU with coordinates i = 66; j — 41 (cells 
are numbered beginning from the left bottom corner of 
the grid). But evolution of Pmr is qualitatively different 
from the Pml case. At n — 500 (right column, row b at 
Figure 10) Pmr is fairly asymmetric. At n = 1256 one 
can see that Pmr takes form of a RSW with Qt = +4 
(right column, row c at Figure 10). But this RSW is 
unstable — the attractor for Pmr is RSW with Qt = +3 
and period T — 2308 (right column, row d at Figure 10). 

This example shows that hypersensitivity of TLCA to 
changing of initial conditions may lead not only to full 
"breaking" of the slightly imperfect symmetry, but to 
restoring of such imperfect symmetry to the perfect one. 
Besides of C v or Ct cases, a similar competition of sym- 
metry "breaking" and restoring is possible for the cases 
of C 2 , C 4 (both under GBC-1 or GBC-2 conditions). If 
starting patterns consists of multiple symmetric parts, 
then similar phenomena will appear at least under GBC- 
1 conditions (for perfect mosaic pattern and the GBC-1 
conditions, a regular vortex grid is obviosly possible). 

Besides of these broken/restored symmetry phenom- 
ena, there are another group of effects, which are con- 
nected with competition of left/right spinning of RSWs. 



4 This is, of course, not a "Cheshire cat" of C. R. Tompkins (a 
two-level block predecessor, which is known in CGL). 



We already described such phenomena in Subsection llV Bl 
for cases of random starting patterns. One can mistak- 
enly conclude that unpredictability of direction of RSW 
spinning is the consequence of randomness of inputs in 
general. But in the next Subsection we will show that the 
main source of unpredictability of handedness in TLCA is 
hypersensitivity to initial conditions, when minimal pos- 
sible change of starting pattern may reverse direction of 
RSW spinning (both at attractor and at transient stages 
of evolution). 



E. Hypersensitivity to initial condition and 
unpredictability of left- or right-handed forms of 
vorticity 

Changing of a single CAU may lead to strongly ex- 
pressed mutation of spatio-temporal structures in TLCA, 
as it was demonstrated in the previous Subsection. Such 
hypersensitivity is an inherent property of very wide class 
of nonlinear multicomponent systems with growing spa- 
tial disturbances and multiple attractors. Some more or 
less close analogs one can find in such well-known area 
of technology as electric power grids, which already has 
been mentioned within the framework of chimera states 
consideration in Subsection llV CI (on cascading effects in 
power grids and other real- world and articicial multicom- 
ponent systems see also [3(| and references therein). 

Now we will investigate mutations which produce re- 
versing of handedness of attractors reached during TLCA 
evolution in conditions of hypersebsitivity to initial con- 
ditions. Let the configurations of ground-state CAUs in 
starting patterns at Figure 11 (Pg[ at top of left column 
and Pgfl at top of right column) are perfectly symmet- 
ric (in distinction to starting patterns Pf/r, an< l ^mr 
shown at Figure 10). So asymmetry in both P^ L and 

*SR (Pig ure 11) is determined only by slightly differen- 
tiated populations of two non-ground levels. 

Namely, left-column starting pattern Pg^ at Figure 11 
is not fully mirror-symmetric, because ten black cells at 
the top of the excited left part of P^ are changed to 

ten gray cells at the top of the excited right part P^ L - 
At n — 2010 pattern P$l takes form of a RSW with 
Qt = +4 (left column, row c at Figure 11). But this 
RSW is unstable — the attractor for Pmr is RSW with 
Qt = +3 (left column, row d at Figure 11). 

The right-column starting pattern P^ at Figure 11 
differs from the left-column starting pattern only by 
the single CAU with coordinates i = 66; j = 41 (similarly 
to that for Figure 10, but black cell at site (66; 41) is 
changed by gray one). Evolution of Psr. is qualitatively 
different from the Psl case and leads to an attractor 
having another \Qt\ and opposite sgn(QT), namely to 
RSW with Qt = —2 (right column, row d at Figure 11). 

This last result illustrates that hypersensitivity of 
TLCA to changing of initial conditions is a cause of un- 
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predictability of the direction of a RSW spinning even 
for the robust vortices. Minimal possible disturbance of 
the initial conditions (by changhing of the state of a sin- 
gle CAU in a starting pattern) trigger the sign of Qt, 
not to speak of the Qt modulo. In other words, a re- 
sult of self-organization of a dynamically stable world in 
TLCA cannot be forecasted even for almost identical ini- 
tial conditions. Left-handed or right-handed life in the 
framework of the model under consideration is a result 
of pure accident despite of extrmely simple, absolutely 
deterministic underlying algorithm. 

In conclusion of this computer-experiment Section it 
will be marked out the following. Rich phenomenology 
of TLCA evolution, illustrated in this Section by several 
important examples concerning RSWs 5 , demonstrates 
that this kind of deterministic mapping is one of the sim- 
plest, fully discrete, bounded (except the case GBC-3) 
system possessing complexity of very high degree. Brief 
list of observed complex phenomena includes: 

• Transient spatio-temporal chaos (including long-term 
irregular evolution of complicated RSW structures); 

• Competition of left-handed and right-handed RSWs 
with chaotic alterations of Q T nt ° SY ' and unexpected hand- 
edness of the winner RSW; 

• Incongruous or fully antagonistic spatio-temporal 
structures (chimerae); 

• Amazing symmetry-connected properties of TLCA 
caused by hypersensitivity to initial condition, including 
full restoring of (initially broken) symmetry of evolving 

TLCA pattern. Reversing of sgn ^Q T Dtcg1 ^ with chang- 
ing of |Qy ntesr '| (mutations) caused by modification of a 
single CAU state only at the beginning of evolution. 

In the Section we discuss relationship between sim- 
plicity of the TLCA model and complexity of its be- 
haviour as well as some other issues concerning self- 
organization in CA. 



V. DISCUSSION 

A model of an active medium based on cyclic transi- 
tions in a three-level system was introduced for the first 
time by astrophysist A. D. Thackerey in 1930's [78j. In 
essence, this was a model of resonatorless laser system, 
but proposed about 20 years before the beginning of the 
laser era (and about 10 years before work of N. Wiener 
and A. Rosenblueth l90l on excitable systems and work 
of J. von Neumann [§2| on cellular automata). 

Similarity of class-B lasers and excitable (chemical) 
systems was revealed and investigated by C. O. Weiss, 
K. Staliunas and their co-workers [IS El El El with- 



out using of CA or any other discrete mapping of this 
kind. Are algorithmically simple (but having phase space 
of huge dimension) CA models of class-B lasers and/or 
laser-like systems appropriate for at least qualitative 
study of these active devices? Are there advantages of 
CA in modeling of laser (phaser) dynamics at all? 

There are many publications on CA and other algorith- 
mically simple models of complex systems, which may be 
joined under the common title: "Complex worlds from 
simple rules?" (this is the real title of the shortest, 
to our knowledge, paper on this subject — see note of 
U. Yurtsever [93). Criticism of such publications is di- 
rected against simplicity of CAU operation, but it usually 
does not take into account emergent co-operation of large 
or very large quantity of CAUs in an intrinsically unsta- 
ble system with random (or randomly perturbed) initial 
conditions. 

It is obvious that discrete, spatially bounded deter- 
ministic system, e. g. CA (ZM or TLCA in partic- 
ular) may have only two types of attractors, both of 
which are regular (as it was already stressed in Sec- 
tion [Wj. In countrary, continuous deterministic system 
may have a great variety of irregular attractors [70J or 
even much more complicated (turbulent) spatio-temporal 
structures 6 . Algorithmic complexity of models simulat- 
ing continuous systems is usually higher (by many rea- 
sons, including purely technical) . Introducing of stochas- 
tic terms in a (classic) model brings new sources of com- 
plexity. And quantum-mechanical models in most cases 
has the highest level of complexity. This hierarchy of 
complexity is built by the principle of algorithmic com- 
plexity. But algorithmic or Kolmogorov complexity is not 
the unique tool and it possibly is not an adequate tool in 
the area under consideration |72| . So (in countrary to the 
U. Yurtsever [93 point of view) "cellular automata are 
one of the more popular and distinctive classes of models 
of complex systems" (see [T^, page 20). 

Results of the present work demonstrate complex, un- 
predictable behaviour of algorithmically simple system 
containing mesoscopic quantity of discrete elementary 
units (usually 10 4 — 10 6 ). We would like to stress that 
we deal with high-dimensional system, which cannot be 
reduced to any averaged, low-dimensional one. Distinc- 
tions between "standard" system having dimension of 
phase space d p hase = 3 and a system with d p hasc = 
10 4 — 10 6 is, of course, very significant circumstance [B| 
in modeling of active systems containing discrete inter- 
acting particles. Phaser system is one of such systems. 
Real phaser contains macroscopic quantity of particles 
(dphasc > 10 16 -10 17 ) |HHIiiSlBlll,butweho P e 
that mesoscopic system is more appropriate for model- 
ing of spatio-temporal chaotic ( "turbulent" ) motions , 



5 Besides of RSWs, there are several qualitatively another types of 
robust spatio-temporal structures in TLCA, and a lot of fragile 
ones. They will be discussed in separate publication. 



6 Some intermediate cases are represented by discrete, spatially 
unbounded systems (e.g. CA of GBC-3 type, see Section ITTT1. 
where true irregular spatio-temporal structures may appear. 
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self-organization and slow self-detunings 0> 0, 03, , 
inversion collapse pM l43l | and other nonlinear phenom- 
ena in this multi-particle dissipative system 7 . 

In particular, transient spatio-temporal chaos observed 
in TLCA may be collated with complex, "turbulent" be- 
haviour of autonomous phaser generator 01 . On the 
other hand, experimentally revealed coexistence of reg- 
ular and irregular motions in ruby phaser [4fij may be 
interpreted as chimera states 0, observed in TLCA 
as transient phenomenon (see Figure 9). 

Generally, TLCA model describes oscillatory collective 
states in system of non-oscillatory CAUs, including local 
self-synchronization of oscillations caused by diffusion of 
excitations, and global synchronization (or, alternatively, 
full collapse of excitations) as final state of the system 
evolution. In the first case a robust RSW or more com- 
plicated but strictly periodic in time structure is formed 
being a cyclic attractor of TLCA at several initial con- 
ditions — see e. g. Figure 3. This is an analog of limit 
cycle, i. e. regular attractor known in lumped dynamical 
systems. In the second case, when all non-ground CAU 
break out at finish, the full freezing of TLCA takes place 
- see e. g. Figure 4. This is an analog of point-like 
attractor in lumped dynamical systems. Both attractors 
may coexist for the same set of CP, as it is for the cases at 
Figure 3 and Figure 4. Such coexistence of final states is 
not the full analog of generalized multistability, which is 
well known for lumped dynamical systems, but branching 
of self-organized states in TLCA is evident. 

It is interesting to discuss shortly the question of degree 
of self-organization in the systems of CA type [23 . It is 
obvious intuitively that final robust RSW in a TLCA 
(Figure 3) is a result of high-degree self-organization in 
TLCA. But it is not obvious intuitively that absolutely 
frozen final state (Figure 4) of TLCA is self-organized. 
At the same time, both routes to final states at these 
figures are very similar. This or that attractor is reached 
due to predefined (but deeply hidden) route dependent 
on combination of TLCA rules, control parameters and 
initial conditions. Slight changing of initial conditions 
frequently leads to another attractor, but the evolution 
itself may by qualitatively the same. If evolution is highly 
bottlenecked then attractor may not be reached at all 
during the whole time of numerical or real experiment. 

From this point of view a final state itself may be of 
minor interest, more important is a transient dynamics. 
The work of J. P. Crutchfield and K. Kaneko |18j . which 
was already cited in Section llVl was titled: "Are Attrac- 
tors Relevant to Turbulence?" . We may slightly reformu- 
late this question in the context of our study: "Are At- 
tractors Relevant to Transient Spatio- Temporal Chaos?" . 
The answer is "Yes" if an attractor may be reached for 
a reasonable time t at tr (in fully discrete and bounded 



7 Phenomena observed in driven phasers |2.1 143t I4SL l4fll I "ill 

need modification of our TLCA model to a non- autonomous one. 



system i a tt r is always limited by quantity of all possible 
states of the system). But the answer is "No" if t at ti 
exceeds any possible duration of an experiment. In this 
case a system with transient spatio-temporal chaos can- 
not be distinguished from true chaotic system without 
additional testing. 

As a matter of fact, there are some intermediate classes 
of phenomena "at the edge between order and chaos" 
which may appear in bounded discrete deterministic sys- 
tem with large phase space. And self-organization sce- 
nario which includes super-slow, bottlenecked, chaotic- 
like stages is a signature of dominance of such an in- 
termediate class of system dynamics in numerical exper- 
iments. It is important to define not only qualitative 
criteria of self-organization, but quantitative ones too. 
Really, having limited time and computer capacity, one 
cannot reach final self-organized state for a system with 
huge dimension of phase space. Irreducibility of CA or 
another discrete mappings does not permit direct fore- 
casting of the system future without direct computation. 
So the computed part of transient process is the single 
source of available information of our fully deterministic 
but partially determined system. 

Criteria of self-organization introduced during last 
decades are of great interest. The brightest of them, 
absolutely counter-intuitive example is suggestion of 
Yu. L. Klimontovich (3?|], which may be formulated 
shortly as follows: Turbulence is more organized state 
than laminar motion. This suggestion is based on us- 
ing of normalized entropy of nonequilibrium state. Para- 
doxical approach of Klimontovich is actually deep and 
have definite perspective for development in nonequilib- 
rium thermodynamics. But it is inadequate for very wide 
class of (self-)organized complex system including high- 
dimensional system of CA type. 

Nonequilibrium thermodynamics is very capricious 
thing, it may give both excellent and unappropriate re- 
sults for the same physical object under different con- 
ditions. Our previous experiments on phaser systems 
confirmed the need of careful analysis of applicability 
of thermodynamical approach to concrete cases. Very 
good qualitative and even quantitative interpretation of 
experiments on phaser amplification is contrasted to the 
case of phaser generation, where spin-temperature model 
sometimes gives very bad description of phenomena ob- 
served (and such models are failed to produce any heuris- 
tics of phaser generation). In an amplifier, there are no 
self-organized states at all, and saturated populations of 
levels admit avera ged description (by spin temperature 
in rotating frame) |24l — in contrary to some regimes of 
phaser generator [44j. Using of exact high-dimensional 
CA models is, to our mind, more appropriate for investi- 
gation of self-organized multiparticle systems. 

Recently some interesting results were reported on 
quantitative investigation of self-organization in cyclic 
CA [73, which are a kind of multilevel CA H with SSR 
(see Section Cyclic CA are more expedient for chem- 
ical systems, which haven't real discrete levels. Phasers 
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(as far other active systems based on inversion of energy- 
level populations), in contrary to chemical excitable sys- 
tems, have well resolved, discrete levels (usually three or 
four), so TLCA is an adequate model from this point 
of view. Moreover, our TLCA is of MSR type, so it is 
a flexible tool taking into account relaxation times for 
nonground levels. Application of the approach, proposed 
in [73j, to the TLCA may give a useful tool for direct 
quantitative study of self-organization in this cellular au- 
tomaton and, possibly, for understanding of nonlinear 
phenomena observed in real physical experiments with 
Gigahertz phaser generator [43L 13 I52I IH^ | and similar 
active systems with negligible level of quantum noise. 

And a bit about applicability of the TLCA to modeling 
of optical-range phonon lasers must be explained. Several 
successful attempts of laser-like generation of phonons at 
low-frequency side of Terahertz range were undertaken 
at the end of 1970-th |58|, and now the renewed interest 
to this frequency range is coming into sight |7jj. There 
are two main differences between 10-Gigahertz (usual mi- 
crowaves) and, say, 100- Terahertz (or 1000-Terahertz) 
ranges: (i) about 12 (up to 15) orders difference in level of 
quantum noises, already disscussed in Section [i] and (ii) 
sufficient distinction in distribution of level-population, 
which will be considered here for the case of three-level 
system. 

In an excitable medium of a Gigahertz quantum de- 
vice, all three working energy levels (ground, excited and 
refractory levels in context of the present work) are pop- 
ulated even without inversion and even at liquid helium 
temperatures. Short pulse inversion almost instantly 
changes relative distribution of active centers with dif- 
ferent states, but after the inversion pulse ending, this or 
that part of active centers remains in the refractory state. 
Examples of such distributions are shown at Figure Al. 
In this case the TLCA algorithm is directly applicable 
for modeling of spatio-temporal dynamics in active (ex- 
citable) medium. 

On the other hand, at optical-range frequencies only 
the ground level L\ is populated in equilibrium system. 
Pulse inversion primary populates the highest level Lm, 
but there still are no active centers at intermediate (re- 
fractory) level L\\ during at least first stage of inversion. 
To model dynamics of such the medium correctly, the 
TLCA must be complemented by appropriate branch(es) 
of evolution operator describing pointed critical stage of 
population distribution. This is a future direction of 
modeling of laser-like active systems by CA. 



VI. CONCLUSIONS 

In this work, we fulfill computer modeling of spatio- 
temporal dynamics in large (up to 10 6 particles) dicrete 
systems of three-level CAUs interacting by short-range 
diffusion through one or two channels. Long-time evolu- 
tion (10 5 — 10 7 iterations) of TLCA demonstrates some 
unusual phenomena including highly bottlenecked col- 



lective relaxation of excitations by multiple RSWs with 
variable Q T ; competition of left- and right-handed RSWs 
with alteration of integral Qt; long lived chimera states, 
i. e. coexistence of regular and chaotic domains in evolv- 
ing patterns of CAUs; branching of TLCA states with dif- 
ferent symmetry which, in particular, leads to full restor- 
ing of symmetry of imperfect starting pattern. Most of 
observed processes of TLCA evolution may be attributed 
to transient high-dimensional spatio-temporal chaos in a 
deterministic excitable system of three-level active units. 
Super-slow collective motions, spatio-temporal transient 
chaos and coexistence of regular and chaotic structures 
in TLCA qualitatively describe part of real experimental 
data on phaser (microwave phonon laser) dynamics. 



APPENDIX A: ENERGY LEVELS AND 
MICROWAVE-FREQUENCY TRANSITIONS IN 
Ni 2+ : A1 2 3 SPIN SYSTEM 

Free 3d 8 -ion Ni 2+ has ground term 3 F (see, e. g. pj 
S HI]). This term consists of (2£ + 1) x (25 + 1) = 21 
levels, which are unsplitted for a Ni 2+ ion in free space 
(Figure A2). Here £ and S are the orbital and the spin 
quantum numbers respectively: C = 3, S = 1. 

In corundum (AI2O3) crystal matrix doped by nickel, 
each Ni 2+ ion is under action of static six-coordinated 
trigonal electric field (called crystal field 0,E|) with the 
highest-symmetry axis O3. This third-order axis coin- 
cides with the optical axis of corundum crystal matrix. 
In the Appendices EI El we will use the Cartesian coor- 
dinate system with the applicata O z directed along the 
axis O3. 

The electric field of the trigonal symmetry splits the 
ground term 3 F into orbital singlet 3 A2 and two orbital 
triplets 3 Ti and 3 T 2 (Figure A2). The lowest 3 A 2 or- 
bital singlet in non-relativistic approximation has three 
unsplitted spin levels (the orbital singlet 3 A2 is simulta- 
neously the spin triplet), because crystal (electric) field 
in this approximation does not interact with spin (mag- 
netic) moment. 

But this degeneracy is partially takes off by the rela- 
tivistic spin-orbit interaction, described by operator \CS 
P, |ij , where A is the constant of spin-orbital interaction; 
C and S are the operators of the orbital and the spin 
momenta respectively. As the result, spin triplet 3 A 2 
becomes splittcd into low-lying spin doublet and the ex- 
cited spin singlet (Figure A2). The spin singlet for the 
Ni 2+ : AI2O3 spin-system lies at V /h « 39.8 GHz over 
the doublet P, |g , where T>q is so called zero field (more 
rigorously, zero-magnetic-field) splitting of quantum en- 
ergy levels. 

This splitting in Ni 2+ : AI2O3 is about fourfold greater 
than T>q in ruby (Cr 3+ : AI2O3) due to sufficient dif- 
ference in |A| for Ni 2+ and Cr 3+ ions: JA(Ni 2+ )| = 
335 cm -1 versus |A(Cr 3+ )| = 87 cm" 1 6]. Accord- 
ingly, spin-phonon interaction in Ni 2+ : AI2O3 is much 
stronger (and longitudinal paramagnetic relaxation is 
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much faster) than in ruby — compare experimental data 
on tensor of spin-phonon interaction for Ni 2+ : AI2O3 pEj 
and Cr 3+ : AI2O3 [13 )• We want to underline close inter- 
connection between static spin-lattice interaction (which 
defines zero-magnetic-field structure of the low-lying en- 
ergy levels) , dynamical spin-phonon processes (which re- 
spond for the resonant interaction of paramagnetic ions 
with microwave ultrasound) and longitudinal relaxation 
in such the diluted paramagnetics. The features of longi- 
tudinal relaxation of Nickel ions in AI2O3 will be consider 
in the Appendix [BJ Now we will continue description of 
the energy levels and microwave-frequency transitions in 
the Ni 2+ : AI2O3 paramagnetic system. 

In nonzero static magnetic field H, which is applied 
along O3, the doublet splits into pair of Zeeman lev- 
els with M = ±1, and the three-level spin system 
Ni 2+ : AI2O3 becomes fully splitted (Figure A2). The 
static magnetic field does not affect the level with M = 0. 
Here {M} = { — 1, 0, +1} is the set of the eigenvalues of 
the scalar operator S z (this operator is the projection of 
the vectorial spin operator S onto applicata): 

§ z \Vm) = M\V M ), (Al) 

where \^m) are the wave functions (eigenfunctions of 
S z ). 

The energy levels Eim) of the Ni 2+ : AI2O3 spin system 
for the case under consideration (H || C 3 ) are as follows 

E {0) = (+2/3)P ; (A2) 
E (+1) = (-1/3)D„ + H\\PbH; (A3) 
E ( _ 1} = (-l/3)D -fig||/Jnff, (A4) 

where g|| is the appropriate component of the effective 
g-factor; (3b is the Bohr magneton; H = |H|. 

From the practical point of view (e. g. for using 
the Ni 2+ : AI2O3 as phaser active system [2j| ^3|) the 
most interest is the case of small magnetic fields: H < 
T>o/H , yu. In this case the following correspondence be- 
tween Ni 2+ : A1 2 3 levels E {M) and TLCA levels E K 
( Section ITTTfl is suggested: 

E {0) Lur, (A5) 
E {+1) & L„; (A6) 
£(-1) & (A7) 

Possibility of such modeling of the Ni 2+ : AI2O3 active 
spin system by the 2C cellular automaton (Section 
is sustained by the obvious selection rules for transitions 
between E^m)'- 

, AM , f 1. fOT (£ ( _ x) - E m ) A (E {0) -> 
[2, for (E {+1) 

The main mechanism of diffusion of spin excitations in 
paramagnetic systems of the Ni 2+ : AI2O3 type is the 



magnetic dipole interaction between active units at res- 
onance frequencies 0, 0, E3- These interactions are 
permitted for \AM\ = 1, and they are forbidden for 
AM = 2. So, there are only two channels for diffu- 
sion of spin excitations in physical system Ni 2+ : AI2O3 
at H || O3 — this is the case of the diffusion of exci- 
tations in the TLCA model fSection 1111(1 . As the result 
we have such the correspondence between £Vm) — > E(m') 
and Lk — ► Lri 

(E {0) - E {+1) ) (Lm - La) ; (A9) 
(E {+1) -» «*• (L u -> ii) ; (A10) 

(£?(_!) -»^(o)) (ii ->£m). (AH) 

For H j| O3 this close correspondence is, generally speak- 
ing, destroyed. The larger is angle a between H and 
03, the more mixing of wave functions ^m} for each 
level takes place. But the case of H || C 3 (in practice 
I a I ;$ 5°) is the single configuration where the phaser gen- 
eration in the Ni 2+ : AI2O3 was experimentally realized 
|43j. An attempt of P. D. Peterson and E. H. Jacob- 
sen [(35j to excite phaser generation in Ni 2+ : AI2O3 at 
a = 76° was unsuccessful despite of large amplification 
of injected microwave ultrasound (some possible reasons 
of this "silence" of an active system Ni 2+ : AI2O3 in ex- 
periments of P. D. Peterson and E. H. Jacobsen [65| were 
discussed in 0|). ^ n this work we assume H || O3, and 
the Ni 2+ : AI2O3 system is modeled by TLCA on the ba- 
sis of (HJl-lUnj 

APPENDIX B: TIMES OF LONGITUDINAL 
PARAMAGNETIC RELAXATION IN Ni 2+ : A1 2 3 
SPIN SYSTEM 

Most of our early microwave experiments on spin- 
phonon interaction and phaser generation in active para- 
magnetic medium Ni 2+ : AI2O3 were fulfilled j43| at 

H <£V /h g[l f3 B , (Bl) 

i.e. under conditions: 

^(+1,-1) « {F {0 ,-i),F [0!+1) ,V /27Th), (B2) 

where F(m,a/') ar e the frequencies of quantum transitions 
in our three-level spin-system: 

F(m,m') = [E ( M) ~ E (M ,))/2nh (B3) 

Typical values of CP in experiments [2^, ^| were as 
follows: H w 0.5 kOe; F {+1 _ x) = 3.0 GHz, F (0 ,-i) « 
41.3 GHz, F ( o !+ i) » 38.3 GHz. At first blush, such the 
system may be described by TLCA with r e <C r r , be- 
cause time of direct one-phonon longitudinal relaxation 
ji(direct) ^ ^ e two-level approximation for each the tran- 
sition) depends on -FW,M') by such the way 

^(direct) k _3 tanh (xF {MiMI) /6) . (B4) 
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Here 9 is temperature of the crystal, x is normalizing 
factor. We use the standard denotation Ti for the lon- 
gitudinal relaxation time in physical systems instead of 
r, which is used for relaxation time in the TLCA model 
system. In this case r r in the TLCA model (Section IIII|) 
is an analog of Ti(Et±x)), and r r is an analog of T±(E^) 
(this analogy is, of course, only qualitative, but it is clear 
from the physical point of view and expedient for numer- 
ical modeling of active systems) . 

For the case T(m,M') = ^(+1.-1) = 3 GHz the inequal- 
ity (>cF(m,m>)/9) 1 takes place at liquid helium tem- 
peratures (8 < 4.2 K). So the time of longitudinal relax- 
ation y.( dlI0Ct ) f or the level (which is the analog of 
the refractory level L\\ in TLCA — see Appendix lAl at 
these temperatures must have such the form: 

T^ ct \E {+1) ) oc F [+ l_ 1} e-\ (B5) 

This two-level approach (which is often used in modeling 
of passive nonlinear systems) leads to wrong conclusion 
that relaxation of the level -E(o) must be much faster than 
relaxation of the level and that the correspondent 

TLCA with r e < T r is a good model for the Ni 2+ : A1 2 3 
physical system 8 . 

In a three-level (or, generally, multilevel) spin system, 
especially with large difference in splittings, there usu- 
ally are several different mechanisms of longitudinal re- 
laxation (direct, Orbach, Raman etc. Q,|(|): 

1/Ti = (1/Tf rect )) + 

+ ^/T^ 0rbach ^) + f]yy>.( Raman h -)-... 

Measurements in Ni 2+ : AI2O3 at low temperatures pM 
Ei l showed that the dominant mechanism of longitudinal 
relaxation for the spin doublet M = ±1 at Fr+x,-i) <C X>o 
is the Orbach two-phonon process Q, ||| : 

r 1 (E (+1) )«Tf rbach )(^ (+1) ), (B7) 

with: 

T (Orbac h ) {E(+i)) ^ W S exp(xPo/0) w 

w A: F ( - 3 +1) exp(xF (0 , + i)/fl) « 7i(£ (0) ), 

where fco is normalizing factor. As the result, mentioned 
early inequality r e <C T r is generally incorrect for mod- 
eling of dynamics of the Ni 2+ : AI2O3 active system (at 
least for H || O3, 9 — 1.8-4.2 K, where the experiments 



Relaxation times r e and r r are, of course, not the precise analogs 
of the correspondent physical relaxation times Ti(_B(q)) and 
Ti (_£?(_)_!)) (the last are not additive etc.). Nevertheless, very 
close correspondence is obviously present here, and not only 
qualitative but quantitative interrelations between physical re- 
laxation times may be used for defining of r e and r r in a TLCA 
model of physical (particularly, phaser) media. 



on phaser generation [25t Elij were fulfilled) . In the case 
under consideration, an approximate equality r e ss r r 
takes place instead of the above strong inequality, as it 
follows from Eqns. (|B7I) and l]B8fl. 

Additional changing of physical relaxation times at 
spin transitions of Ni 2+ : AI2O3 is caused by interaction 
(cross-relaxation and some other mechanisms) of Ni 2+ 
ions with the Jahn- Teller 3d 7 ions Ni 3+ |43j in nickel- 
doped corundum. In strong six-coordinated crystal field 
the single unpaired 3d electron of Ni 3+ ion is at e-or- 
bital, i. e. the ground state of Ni 3+ in corundum is 2 E 
(orbital doublet and spin doublet simultaneously) [74| . 
Non-zero ground state orbital momentum of the Ni 3+ ion 
in the Ni 3+ : AI2O3 system (in contrast to zero ground- 
state orbital momentum for Ni 2+ : AI2O3) is the cause of 
strong electron-phonon interaction and fast longitudinal 
relaxation of the Ni 3+ [74[ . 

The last type of ions always is present in nickel doped 
corundum, because ion Ni 3+ has the same charge as Al 3+ , 
and after usual crystal growing the main part of nickel 
ions is in trivalent state. Reducing of Ni 3+ to Ni 2+ needs 
special technological operations, and such a reducing can 
not be full in principle. 

Interaction of Ni 2+ with Ni 3+ in corundum strongly 
depends on magnetic resonance frequencies F^ M ^ M i^{H), 
this interaction is very sensitive to concentrations of both 
ions in corundum, etc. As the result, the ratios of relax- 
ation times in Ni 2+ : AI2O3 spin-system varies in very 
wide range, depending not only on experimental condi- 
tions, but on a crystal growing technology too. Due to 
these circumstances, the adequate TLCA model of the 
Ni 2+ : AI2O3 active system must have possibilities for 
working with arbitrary ratios r e /r r . 



APPENDIX C: COLLAPSE OF INVERSION 
STATES AND CRITICAL SLOWING-DOWN 
IN A THREE-LEVEL PHASER AMPLIFIER 
WITH BISTABLE PUMPING 

Quantum amplifiers with bistable pumping were pri- 
marily investigated for the simplest case of weak signal 
pH IItI , when amplified field does not affect the active 
medium. The spin system of nonlinear three-level phaser 
amplifier (of a Ni 2 + : A1 2 3 type HHHHSHJ) is sat- 
urated by two microwave fields simultaneously. 

The first of them is the electromagnetic pumping field 
with frequency F pump = i*(o,-i) and normalized ampli- 
tude of magnetic component Y . This microwave mag- 
netic component interacts with the spin system, exciting 
I AM I = 1 transitions between the lower T£(-i) and the 
upper E/m spin levels. 

The second field is microwave ultrasound (called also 
hypersound) with frequency -F s i g nai = ^(+1.-1) and nor- 
malized acoustic intensity J. This acoustic field interacts 
with the spin system, exciting |AM| = 2 transitions be- 
tween the lower E(_u and the intermediate spin 
levels. Such unusual selection rule (|AM| = 2) is caused 
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by the quadratic on S Hamiltonian of spin-phonon inter- 
action [EEJIISII2 f° r iron group ions in crystals. 
Let us normalize the population differences Nim m') = 



N(m') a t the spin transitions £ym) — ► E, 



their equilibrium values N, 



(M,M> 



(m>; 



to 



AT, 



(M,M')|y =0 ,J = 



D 



A< ( 



(-1,0), 



N, 



(+1,-1) 



A<, 



(e) ' 
(-1,0) 



A', 



(e) 

(+1,-1) 



(CI) 



where < D < 1 and — 1 < K < L — 1 (note sign "minus" 
in the definition of K, Eq. IjCll) ). Here K is the inversion 
ratio at the — * -E(-i) spin transition, and the pa- 
rameter L characterizes properties of an active medium. 
Value L depends on the ratio of pump and signal frequen- 
cies -Fpump/^signai, on the times of longitudinal relaxation 
etc. The third population difference A^ +1 ) is univocally 

defined by D and K, because J2m^(M)/^(m)) = 

For linear autonomous quantum amplifiers with weak 
signal (J <C 1) standing-wave pumping becomes bistable 
pHE?! if parameter of cooperativity C is greater than its 
critical value C cr (where C cr is codimension-2 bifurcation 
point), and if Yj, < V < Yj (Yj, and Yy are codimension-1 
bifurcation points). For nonlinear autonomous quantum 
amplifiers bifurcation values C cr , Y^ and Y^ are renor- 
malized by an intense running-wave signal with constant 
amplitude [50l l5l| . and the collapse of inversion state 
becomes possible. 

Phaser amplifier with modulated microwave acoustic 
signal is nonautonomous (at least one CP in our dy- 
namic system becomes time dependent). But the equa- 
tions of motions for nonautonomous dynamic system may 
be transformed to e quiv alent autonomous form in an ex- 
tended phase space |64|. The only requirement for this 
transformation is regularity of modulation. In an ex- 
tended phase space (D,K,£) c M 3 (where ( = ui m t, 
Lo m is the modulation frequency, and t is usual, i. e. 
non-discretized time) the equations of motion for nonau- 
tonomous phaser amplifier with bistable pump |50l l5l| 
transforms to such the autonomous form [54| : 



(C2) 



{T 1 dA/dt = F(AX/u m ,® A ); A(0) = A ; 
\d(/dt = uj m ; C(0) = w m to 

Here A = (D, K); 0,4 is time-independent vector of CP; 



F = (F D ,F K ). We suppose T x = T X {E^) « Ti(£7 (+ i)) 
(i. e. r e k r r ) by arguments of Appendix ITU and the 
inequality characterizing class-B active systems T\ » 
Tficid ^> T% takes place by definition [s^] (see also Section 
Uof the present work). Components of vector F are as 
follows: 



Ed 



1 - D 



1 - K 



Y 2 D 



(1 + 2CD) 2 4L' 
ly2 ° JK, 



#; (03) 



(1 + 2CD) 2 



(C4) 



where J = (l + i; m j sin () Jo; k m j is the coefficient of 
modulation. 

Eqs. l|C2p to (|C4p describe bistability, collapse of in- 
version states, critical slowing-down and other nonlinear 
processes in nonautonomous phaser amplifier (in frame- 
work of lumped, i.e. point-like three-level model). An 
example of critical slowing-down during collapse of in- 
version states is shown at Figure A3. 



APPENDIX D: LIST OF ABBREVIATIONS 

ID, 2D, ... — One-Dimensional, Two-Dimensional, 

1C, 2C, ... — One-Channel, Two-Channel, ... 

BZ — Belousov-Zhabotinskii 

BR — Bogach-Reshodko 

CA — Cellular Automaton 

CAU — Cellular Automaton Unit 

CP — Control Parameters 

CGL — Conway's Game of Life 

CML — Coupled Map Lattice 

GBC — Geometry and Boundary Conditions 

GH — Greenberg-Hastings 

GWR — Generalized Wiener-Rosenblueth (model) 

MSR — Multi-Step Relaxation 

OK2 — Oono-Kohmoto 2D (model) 

RtC — Rays through Clouds (automaton) 

RSW — Rotating Spiral Wave 

SSR — Single-Step Relaxation 

TLCA — Three-Level Cellular Automaton 

WR — Wiener-Rosenblueth 

ZM — Zykov-Mikhailov 
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Figure 1: Typical stages of TLCA evolution. CP set is as follows: r e = r r = h = 50; g = 1; / = 8. Starting 
pattern p(°) = has 100 x 100 = 10 4 CAUs. Black, gray and white pixels denote excited, refractory and ground- 
state CAUs respectively. Numbers under patterns are equal to the steps n of evolution. At n = 3.66 • 10 5 the system 
is already at periodic attractor. 

Figure 2: Divergence of spatio-temporal structures for / = 8 (upper row) and f — 9 (lower row) during 

TLCA evolution. Values of r e , r r , h, g are the same as in Figure 1, but for another starting pattern P(°) = p[? (see 

Figure Al for Pq )■ At n = 3 • 10 5 both systems are at their attractors. Only excited CAUs are shown at this figure 
(black pixels). 

Figure 3: Slow evolution of TLCA and metamorphoses of RSW topology. Starting pattern P(°) = p£?\ 
Values r e ,T r ,h,g are the same as in Figure 1 and Figure 2, but / = 4. Only excited CAUs are shown at this figure 
(black pixels). The system reaches an attractor with Qt = —3 at n ~ 4 • 10 6 . 

Figure 4: Collapse of excitations for TLCA by slow evolution. Values r e ,r r , h, g, f are the same as in 
Figure 3, but starting pattern is another: P^ ^ = P]^ (see the P]^ at Figure Al). Transient time nc — 3385479 is 
the lifetime of nontrivial state of the whole cellular automaton with the pointed starting pattern. 

Figure 5: Developed transient spatio-temporal chaos in TLCA. Set of CP is as follows: r e = 9,r r = 12, h = 
25; g = 1; / = 9 (ZM model). Starting pattern p(°) = P$\ Lifetime n c = 2427063. Black, gray and white pixels 
denote excited, refractory and ground-state CAUs respectively. 

Figure 6: Pre-finish stage of transient process in TLCA with / = 3. All RSWs have the same (by 
magnitude and sign) effective topological charge. Black, gray and white pixels denote excited, refractory and 
ground-state CAUs respectively. The perfect left-handed multi-RSW pattern is shown for n = 1.6 • 10 5 . The CP set 
is as follows: r e = r r = h = 50; g = 1; / = 3. This pattern is unstable. Multiple RSWs with Qt — —1 compete 
strongly after the left-handed pattern is formed, and "the winner takes all" at n > 3 • 10 — the single RSW with the 
same Qt — — 1 occupies the whole active medium (not shown). But winner must die — excitation collapses fully at 
nc = 396982 due to collision of the winner's core with boundary (not shown). 

Figure 7: Reversing of sign of effective topological charge during TLCA evolution. Part 1: Almost 
purely right-handed pattern. Step n — 1.48T0 6 of evolution of starting pattern P(°) = Pg C with Mx = My = 900 

(pattern P^ itself is not shown here). The CP set is the same as at Figure 6. Only excited CAUs are shown by black 
pixels. Multiple RSW with Qt — +1 forms right-handed pattern. One may expect that winner will be right-handed 
too. But the subsequent evolution of the system is very surprising, see Figure 8. 

Figure 8: Reversing of sign of effective topological charge during TLCA evolution. Part 2: The single 
left-handed RSW is the result of evolution of the right-handed pattern of RSWs. This is one of the 
snapshots (namely at step n = 2 • 10 6 ) of subsequent evolution of the pattern shown at Figure 7. The winner RSW, 
shown here, has Qt = —1, i- e. the system becomes left-handed despite of their long right-handed previous life. During 
further evolution (n > 2 • 10 6 , not shown here), the winner's core is moving irregularly across the grid, it collides with 
the grid boundary, and full collapse of excitations takes place at nc = 2260964 with final stage by scenario shown at 
Figure 4. 

Figure 9: Transient chimera states, i. e. coexistence of periodic and aperiodic spatio-temporal struc- 
tures in TLCA (demonstrated by sequences of generalized Poincare cross-sections). Starting pattern 
p(0) = p(o) Excite d CAUs are black (TRUE), and non-excited CAUs (both refractive and ground-state ones) 
are white (FALSE). Parameters r ei T ri h,g are the same as in Figure 1, but / = 5. Central column of the Figure 
- usual sequence of TLCA patterns, left column — forward-updated sequence of ascending Poincare cross-sections 
V^ + \ni, An, k p ), right column — backward- updated sequence of descending Poincare cross-sections V^(n2, An, m q ). 
Here n\ — 154000; n-i ~ 338000; An = 1000; fc max = ni max — 184 (i. e. ending point for the whole sequence of cross- 
sections is the starting point for pV) an d vice versa). The TLCA at n 2 = 338000 is already at the attractor. 
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Figure 10: Symmetry restoring (left column) and "breaking" (right column) in TLCA. Black, gray and 
white cells denote excited, refractory and ground-state CAUs respectively. Dimensions of the patterns are Mx = 75; 
My = 50 (grid lines are shown for clarity). The CP set is as follows: r e = 5; r r = 7; g = 0.3; h = 3; / = 9. 
Denotations of the rows corresponds to the same steps of evolution for the left and the right starting patterns: 
a — » n = 0; b — > n = 500; c — > n = 1256; d — > n = 10129. Starting patterns differ by the single CAU with i = 66; 
j = 41 (cells are numbered beginning from the left bottom corner of the grid). This CAU in the left starting pattern 

PjSx is at level Lux (excited state, black cell). Analogous CAU in the right starting pattern P^r i s a * level L\ 
(ground state, white cell). 

Figure 11: Hypersensitivity to initial condition as cause of unpredictability of left- or right-handed 
vorticity. Dimensions of the patterns and the CP set are the same as at Figure 10, but starting patterns Pg°£ (left 

column, row a) and Pg^ (right column, row a) are slightly changed comparatively to P^f L and P^L- Denotations 
of the rows corresponds to such the steps of evolution: a ^ n = 0; b ^ n = 500; c — * n = 2010; d — > n = 10118. 
Starting patterns Pg2 an d Pg» differ by the single CAU with the same coordinates i = 66; j — 41 as for Figure 10. 
This CAU in the left starting pattern P^2 is at level Lm (excited state, black cell). But analogous CAU in the right 
starting pattern Pg ^ is at level L\\ (ground state, gray cell). 

Figure Al: Samples of starting patterns. The population of excited CAUs (relative quantity of black pixels) is 
equal to 0.4 for all these patterns. This is slightly less than population of ground-state CAUs (white pixels), but more 
than twice greater than population of refractory CAUs (gray pixels). Resulting inversion of populations at Lm <-> L\\ 
transition corresponds to an alternative variant of the phaser medium activation scheme (in the original scheme of 
the Ni 2+ : AI2O3 phaser (2^,01, the analog of L\\ <-> L\ transition was inverted). 

Figure A2: Energy levels of divalent nickel in corundum. This is the simplest system by which microwave 
phonon laser (phaser) amplification, generation and inertial self-focusing were realized at Gigahertz-range frequencies. 
At non-zero static magnetic field, a collection of such active centers (interacting by d — d mechanism) is a kind of 
excitable three- level system with two-channel diffusion of excitations. See Appendix for explanation of the energy 
levels structure and underlying mechanisms of the levels splitting. For relaxation properties of this system at low 
temperatures see Appendix IB1 

Figure A3: Critical slowing down and collapse of inversion states in a lumped three-level class-B active 
system. Various-time oscillating routes to collapse are shown for the inversion ratio K, see Eqns. I|C2|1 to l|C4|l in 
Appendix |0 The invariant subset of CP (i. e. unchanged components of &a) is as follows: C = 5; L — 2; Y = 6; 
Jo = 0.15; k m j — 0.7. The control parameter uj m is varied relatively to relaxation time T\. Parameter "n" enumerates 
dependencies K(t) for different values of w m Ti, namely: (n = 1) -> {u) m T\ = 0.1); 2 -> 0.15; 3 -> 0.19; 4 -> 0.199; 
5 -> 0.2; 6 -> 0.2004; 7 -> 0.2005; 8 -> 0.20053. Initial conditions: D(t = 0) = 0.01; K(t = 0) = 0.7; C(t = 0) = 0. 
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